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Abstract 

Photon trajectories in incoherent radiation trapping for Doppler, Lorentz and Voigt hne shapes 
under complete frequency redistribution are shown to be Levy flights. The jump length (r) dis- 
tributions display characteristic long tails. For the Lorentz line shape, the asymptotic form is a 
strict power-law r~^/^ while for Doppler the asymptotic is (Inr)^^^^. For the Voigt profile, the 
asymptotic form has always a Lorentz character, but the trajectory is a self-afhne fractal with two 
characteristic Hausdorff scaling exponents. 

PACS numbers: 02.50.Ey, 32.80.-t, 32.70.-n 
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Interest in distributions with long tails has increased over the last years, with intensive 
search for such laws in real physical systems. Most of the basic work on these distributions 
was carried out by Levy in the 1930s but it was only relatively recently that these dis- 
tributions were shown to be aDplicable to the description of a number of physical, biological 
and social phenomena B S, fl 0, including particle motion in turbulent media ^M, 
anomalous diffusion in microheterogeneous systems I?)], chaotic transport in laminar fluid 
flow 8], the albatross flight and frequency fluctuations of chromophores isolated in glassy 
environments jl^. The purpose of this Letter is to show that photon trajectories in in- 
coherent radiation trapping, a basic and common phenomenon in atomic and atmospheric 
physics, and in astrophysics, flt in the category of Levy (superdiffusive) flights. 

Although "the general trend nowadays is to put Levy-ty pe a nomalous (super) diffusion on 
a similar footing with normal, Brownian-type, diffusion" IjJ, the actual physical systems 
invoked to justify the practical importance of Levy flights in physical space tend to be 
somewhat marginal in physical phenomena 0|. This contribution aims to prove that the 
well known case of classical incoherent radiation trapping is one of the simplest and best 
characterized Levy flight found up to now. In addition, radiation trapping has a major 
role in fluorescent lamps and this could make radiation trapping the economically most 
relevant concretization of a Levy flight. In retrospect, the 1990 Ott et al. contribution 
on the superdiffusion in microheterogeneous systems can no longer be considered the flrst 
experimental realization of a random flight with inflnite moments, given that it was preceded 
by many studies of radiation trapping. 

Radiation trapping of energy is important in areas as diverse as stellar atmospheres Q| , 
plasmas and atomic vapors luminescence 2, terrestrial atmosphere and ocean optics, 
molecular luminescence infrared radiative transfer and cold atoms Q]. In these op- 
tically thick media, the emitted radiation suffers several reabsorption and reemission events 
before eventually escaping to the exterior; the radiation is said to be imprisoned or trapped. 
Atomic radiation trapping is also known as imprisonment of resonance radiation, line trans- 
fer, radiation diffusion or multiple scattering of resonance radiation. 

The flrst quantitative theory for atomic radiation trapping was presented in the 1920s by 
Compton and Milne, that have developed a modifled diffusion equation for the (frequency) 
coherent spreading of excitation. It was only in 1932 that the frequency redistribution 
between absorption and reemission was taken into account by Kenty He considered 
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a Doppler spectral distribution and arrived to the unexpected result that, for an infinite 
medium, the diffusion coefficient would be infinite. This was the first realization of the 
fundamental fact that all moments of the jump size distribution are infinite. Kenty's re- 
sult shows that a diffusion-type equation is not valid for radiation trapping with frequency 
redistribution effects. Nevertheless it was only in 1947 that Hosltein and Biberman inde- 
pendently proposed a Boltzmann-type integrodifferential equation which remains the 
starting point of the vast majority of radiation trapping models jisl ]. 

Consider the case of inelastic scattering where, as the result of reabsorption-reemission 
events there is a photon frequency redistribution in the lab reference frame. The frequency 
distribution of the emitted photons is given by the emission spectrum O (x) . The absorp- 
tion probability of a photon with frequency x at a given distance from the emission point 
depends on the absorption spectrum $(a;), and is given by p{r\x) = $ (x) e-*(^)^ (Beer- 
Lambert law) where r is the opacity or optical density and is a dimensionless distance. 
Radiation trapping can then be envisaged as a random flight in physical space with spectral 
shape dependent jump size distributions. The jump size distribution takes into account the 
absorption probability for all possible optical emission frequencies, hence 



p{r) = / Q{x)p{r\x)dx (1) 

J — oo 

Eq. dD) fully characterizes the spatial aspects of the random flight. 
The moments of this distribution are 

{r") = n\ / ^^dx (2) 
(x) ^ ^ 



and can be shown to be inflnite for all physical 



reasonable atomic emission and absorption 



spectral distributions as concluded by Holstein 3|. His original analysis only included 
classical incoherent trapping, in which the emitting state is statistically unrelated with the 
absorbing one and there is therefore Complete Frequency Redistribution (CFR). In this case 
the absorption and emission spectra are identical and Eq. ((T)) reduces to 



p [r) 



<l>'(x) e-*(")"dx (3) 



Nevertheless, it can be concluded from Eq. (0) that the case of Partial Frequency Redistri- 
bution (PER) is also characterized by an inflnite moments Levy statistics. The case of PER 
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is especially important in an astrophysical context. Neither case of complete coherent nor 
complete incoherent scattering is achieved exactly in stellar atmospheres and it is then nec- 
essary to consider the photon redistribution and to calculate redistribution functions which 
will give 6 (x) ^12j|. In the usual laboratory conditions, vapor densities are high enough 
for CFR to apply [l3|. This will be the case considered here. In two-level CFR atomic mod- 
els both absorption and emission spectra can be described by Doppler - (x) = ^e~^^ 
Lorentz - = or Voigt - <^v (x) = ^ a^+^I-uf ^ spectral distri- 

butions. X is a normalized difference to the center of line frequency and a is the Voigt 
characteristic width. 

We now consider the asymptotic approximations valid for large jump sizes. A random 
flight in which the probability density of jump lengths is given by 

with /i < 2 is a self-similar random fractal with fractal dimension /i. It is called a Levy 
flight after Mandelbrot [17], and defines a broad distribution for which all the moments of 
order not smaller than /i are divergent. 

If $ (x) is substituted for (x) in Eq. Q above, it is found that (r) = oo, whatever 
the spectral lineshape used, as long as $ (x) is nonzero for large |x|. In this way, /i < 1 for 
any $ (x) . 

In order to find the specific value of /i for the spectral distributions mentioned, we begin 
by rewriting Eq. (jH} as 

P(^) = ^ (5) 

where 

/ + 00 
(1 - e-*(")'-) dx (6) 
-co 

When the line shape is gaussian (Doppler), the integrand in Eq. © approaches a square 
wave form for large r, with inflection points at — xq and xq with $ (xq ) r ^ 1. Hence, xq = 
VTnr for large r, J (r) ~ 2 Xq = 2 y/hir and therefore, from Eq. (0), 
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Assuming an homogeneous scaling law one arrives at an effective /i = 1 + 
1/2 ln{ln{r))/ln{r) which goes to 1 for r oo: although the asymptotic of Eq. (0} is 
only approximately valid, one can nevertheless classify Doppler trapping as a strict Levy 
flight with fi = I with all the moments of the jump distribution being infinite. 

When the line shape is Cauchy-like (Lorentz), the integrand in Eq. (jH)) can be simplified 
to 1 — exp (— $ (x) r) ^ 1 — exp (— for large r and J (r) becomes 2 ^/r. Therefore, 

~ 4^ (8) 

r— >oo 

and the asymptotic distribution is a Levy flight with /i = 1/2. 

The Voigt distribution is asymptotically coincident with the Lorentz distribution, and 
therefore has also yU = 1/2 for any value of its parameter a. 

Eqs. (jZj) and (jH)) were already implied in Holsteins's 1947 contribution (original equations 
are, up to some minor terms, integrated versions of the equations in this Letter), a fact 
hitherto unnoticed in the literature. 

Fig. ^show the jump probability obtained from direct numerical integration of Eq. Q. 
A linear fit gives for the characteristic "tail index" (Hausdorff box counting) fractal dimen- 
sion (yu) of the Doppler distribution fi = 1.07 (jump sizes 10^-10^) and for the Lorentz 
case fi = 0.500 (jump sizes 10^-10*^), in agreement with the previous assertions. Fig. C] 
also show that the continuous transformation from Doppler into Lorentz-like spectra as the 
Voigt a width parameter changes from zero into infinite values does not manifests itself in a 
continuous change of the effective, r dependent, /i values. There is an abrupt change instead 
at a jump distance which scales approximately as 1/a. This is expected as the asymptotic 
expansion is related to the most extreme values of the wings of the spectral distribution 
and these change abruptly from Doppler to Lorentz-type asymptotics for a values as low as 
10~^ (see Fig. 121). The data in Figs. [T] and El allows one to define the Voigt radiation trap- 
ping trajectories as a self-affine fractal with two different scaling exponents which manifest 
themselves at different lengthscales. 

Consider now Fig. El which shows single excitation trajectories in a 3D infinite 
medium (compare Mandelbrot's well known figure Q]). The Lorentz and Doppler cases 
display the two qualitative features characteristic of Levy flights: (i) the longer pathlength 
jumps, although much less common, are of paramount importance to the overall spreading 
of excitation and (ii) self-similar behavior. There is an hierarchy of clusters formed at dif- 
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FIG. 1: Jump size distribution for CFR Doppler, Lorentz and Voigt spectral profiles. From bottom 
to top: Doppler, Voigt with a = 10""^, 10"^, 0.01, 0.05, 0.1, and Lorentz. 




FIG. 2: Doppler, Lorentz and Voigt spectral profiles. 

ferent lengthscales but with similar topology. The set of visited points constitute a fractal 
of characteristic dimension /i for the Doppler and Lorentz distributions . The Voigt 

case shows its self-affine nature. The trajectory topology changes from a Lorentz character 
into a Doppler one for smaller distances. The overall topology of the trajectory as a whole 
is of course dictated by the higher Lorentz scaling. 

Up to now we have considered frequency redistribution at each scattering event. However, 
for high opacity two-level systems there are presumably many elastic scattering events, 
before an inelastic scattering event occurs. It is therefore important to consider at least 
qualitatively the influence of elastic scattering. Elastic scattering events will fold up the 
excitation trajectories of Fig. IHl and the more the higher the ratio of elastic to inelastic 
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scattering probabilities. From this argument alone we will expect that an increase in the 
elastic to inelastic probabilities ratio will lead to an increase in the fractal dimension, as it 
will approach the 3D brownian motion with fractal dimension 2. 

The discussion thus far has been restricted to the case of Levy flights. For the vapors 
studied in the laboratory the time- of- flight of in-transit radiation is negligible compared to 
the waiting time between absorption and reemission, thus rendering the Levy flights formal- 
ism especially adequate. From the set of visited points one can obtain spatial distribution 
functions of the excited density which fully characterize the spatial aspects of radiation mi- 
gration. The temporal evolution can be factorized from the spatial part (Levy flight) and 
then separately handled However, for interstellar gases the time between scattering 

events might be large compared to the absorption-reemission times. Because the speed of 
light is finite, a Levy walk modification of the flights here presented will be more appropriate 
in an astrophysical context jof. 

The theoretical models for incoherent atomic radiation trapping are based either on the 
Holstein-Biberman multi-exponential mode expansion or on the so-called multiple scattering 
representation In the last case, the temporal evolution is known analytically and the 

spatial distributions can be separately obtained, either directly from the Eq. or from 
its asymptotic approximation Eq. (0} 18] . It is within this theoretical framework that the 
results reported in this Letter are most useful. On the other hand, the individual terms in 
the Holstein expansion have no direct relation with the one for n^^) step jump probabili- 
ties. Although both approaches are ultimately equivalen t [l4 , we follow the multiple 
scattering approach, preferable for systems of low opacity 14]. Also, the multiple scattering 
representation has a simple interpretation, since each term corresponds to a specific gener- 
ation of excited atoms or molecules. In a large number of practical situations the opacity 
is not high enough to warrant the exclusive use of Holstein's slowest exponential mode. On 
the other hand, the multi-exponential expansion remains valid, but its fundamental mode 
cannot be identified with Holstein's high opacity result and should be estimated by the sta- 
tionary mode associated with a non-changing spatial distribution function [3 [3, a point 
often misunderstood in the literature. 

Results of Fig. 01 strictly apply to an infinite 3D medium but in a vast number of experi- 
mental situations the system is finite and the trajectory is truncated before the asymptotic 
Levy expansion is able to manifest itself. Mantegna and Stanley j^] introduced in 1994 
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FIG. 3: Single trajectories of 50 000 jumps each for incoherent isotropic CFR radiation migration 
with Lorentz (top row), Doppler (middle) and a = 0.001 Voigt (bottom) profiles in infinite 3D 
medium. Part A show the whole trajectory while B and C show successive details. The three 
trajectories were obtained with the same random number sequence. 

a class of Levy flights, the truncated Levy flight (TLF), in which the largest steps of an 
ordinary Levy flight are eliminated by a sharp cutoff in its power tail. This work allows 
a reassessment of Kenty's pioneering contribution since he was the first to have used a 
truncation procedure: within the framework of the kinetic theory of gases he considered a 



truncated Maxwell distribution of speeds, wit 



1 a maximum speed corresponding to a free 



path equal to the linear size of sample cell |l5|. Although truncation of the jump size ren- 
ders the moments of the distribution finite, the convergence to a Gaussian can be extremely 
slow and the random walk can exhibit anomalous behavior and multi-scaling properties in 
a wide range before convergence 0, 0|- Moreover, trajectory truncation is more complex 
to handle and does not create convergence to the Gaussian in usual situations. 

In this Letter, it was shown that all photon trajectories arising from incoherent two-level 
complete frequency redistribution trapping are super diffusive Levy flights with /i < 1. In 
particular, for the Doppler lineshape = 1, whereas for Lorentz and Voigt fi = 1/2. 
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